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Abstract. We investigate the effect of cooperative interactions on the asymmetric 
exclusion process, which causes the particle velocity to be an increasing function of 
the density. Within a hydrodynamic theory, initial density upsteps and downsteps can 
evolve into: (a) shock waves, (b) continuous compression or rarefaction waves, or (c) 
a mixture of shocks and continuous waves. These unusual phenomena arise because 
of an inflection point in the current versus density relation. This anomaly leads to a 
group velocity that can either be an increasing or a decreasing function of the density 
on either side of the inflection point, a property that underlies these localized wave 
singularities. 
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1. Introduction 

The asymmetric exclusion process (ASEP) [H El El HJ E] represents an idealized 
description of transport in crowded one-dimensional systems, such as traffic |6l [71 E], 
ionic conductors |L9j, and RNA transcription [10]. In the ASEP, each site is either 
vacant or occupied by a single particle that can hop at a fixed rate to a vacant right 
neighbor [Il[2l[3lll]. Although simply defined, this model has rich transport properties: 
for example, density heterogeneities can evolve into rarefaction or shock waves [4j , while 
an open system, with input at one end and output at the other, exhibits a variety of 
phases as a function of the input/output rates p!Tl \T2\ [T3] . 



1 1 X 




Figure 1. Cooperative exclusion. A "pushed" particle (red) — one whose left neighbor 
is occupied — can hop to a vacant right neighbor with rate 1, while an isolated particle 
(blue) hops to a vacancy with rate A. 



A fundamental property of the ASEP is the relation J{p) = p(l — p) between 
the current J and density p. Because each site is occupied by at most one particle, the 
average particle velocity v = J/p is a decreasing function of the density. In this work, we 
investigate a cooperative exclusion (CE) model in which the velocity can increase with 
density. This cooperativity leads to unexpected features in the evolution of initial density 
heterogeneities. Such cooperativity occurs, for example, when ants emit pheromones 
that help guide fellow ants along a trail [H]. Another example are multiple buses that 
follow a fixed route. The leading bus picks up more passengers so that the next bus 
moves faster, which causes clustering of buses during peak travel times [15]. At the 
microscopic level, molecular motors can work together to pull a load that is too large 
for a single motor [16]. Cooperativity has even been proposed as a basis for organic 
superconductors [T7] . 

The notion of cooperative interactions that counterbalance the fundamental 
excluded- volume interaction is implicit in Ref. [7], as well as in [HI [19]. These latter 
publications investigated an exclusion model with a somewhat less stringent excluded- 
volume constraint than in ASEP. This weaker exclusion gives rise to an effective 
cooperativity and thereby to complex density profiles similar to what we find. As 
we shall argue, the existence of these complex profiles does not depend on detailed 
microscopic rules, but is rather a consequence of the underlying cooperative interactions 
between particles. When sufficiently strong, these interactions leads to an infiection 
point in the current-density curve; this feature is the minimum requirement for the 
complex density-profile dynamics. 
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2. Cooperative Exclusion Model 

In the CE model, a particle can hop to its vacant right neighbor at a rate r that depends 
on the occupancy of the previous site (Fig. [1]): 




1 previous site occupied, 
A previous site vacant. 



with < A < 1. When A = 1, the standard ASEP is recovered, while A = corresponds 
to facilitated asymmetric exclusion [20], in which the left neighbor of a particle must 
be occupied for the particle to hop to a vacancy on the right. We pictorially view 
this restriction as a particle requires a "push" from its left neighbor to hop. This 
facilitation causes an unexpected discontinuity in a rarefaction wave in the ASEP pT] . 
More strikingly, we will show that cooperativity leads to shock and rarefaction waves 
that can be continuous, discontinuous, or a mixture of the two. 

These unusual features arise in CE when < A < |, where an inflection point 
in J(p) occurs at p = p/ (Fig. [2]). For p < p/, cooperativity dominates, and J grows 
superlinearly in p. At higher densities, excluded volume interactions dominate, so that J 
grows sublinearly and ultimately decreases to zero. Correspondingly, the group velocity 
changes from an increasing to a decreasing function of density p as p passes through p/. 




Figure 2. Steady-state current as function of density in cooperative exclusion (CE). 
Data are based 10^ realizations with L = 10^ up to t = 10^. The solid curves are given 
by Eq. ([2]). Arrows indicate the locations of the inflection points. 



A configuration of particles on a ring of length L is specified by the occupation 
numbers {ni,...,ni}, subject to conservation J^i'^i = N; here equals 1 if z is 
occupied and equals otherwise. A crucial feature of CE is that the probability for any 
steady-state configuration is a decreasing function of the number k of adjacent vacancies: 
k = Yli=i{^ ~ — Tij+i), with ul+i = rii. To understand how the configurational 
probabilities depend on k, we observe that the hopping of a pushed particle (whose left 
neighbor is occupied) either preserves or decreases the number of adjacent vacancies k 
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(left side of Fig. [3]). Conversely, the hopping of an isolated particle either preserves or 
increases k (right side of Fig. |3]). Since pushed particle hopping events occur at a higher 
rate, configurations with fewer adjacent vacancies are statistically more probable. 



(a) 



(b) 



(c) 



(d) 



Figure 3. [left] Hopping of a pushed (red) particle where the number of vacancy pairs 
is (a) preserved or (b) decreases, [right] Hopping of an isolated (blue) particle where 
the number of vacancy pairs is (c) preserved or (b) increases. 



We now exploit the work of Antal and Schiitz |7] who investigated a dual model in 
which next-nearest neighbor cooperative interactions pull a particle ahead, in distinction 
to the pushing of particles from behind in CE. By the mapping particles O holes, 
the CE and the Antal-Schiitz models give the same probability distribution Pk for a 
configuration with k adjacent vacancies [7]: 

= Z(A) ' 

where Z{X) is a normalization constant. Since A < 1, configurations with fewer adjacent 
vacancies are more probable. Following [7], the steady-state current is 

^l-4(l-A)p(l-p)-l 



J=(l-p) 



2(1 -A)p 



(2) 



in the L — )■ oo limit. The salient feature is that J has an inflection point at a density pi 
for A < I (Fig. [2]). We henceforth restrict our analysis to this domain and determine the 
unusual consequences of this inflection point on the dynamics of initial density steps. 

3. Density Profile Dynamics 

In a hydrodynamic description, the particle satisfies the continuity equation pt + Jx = 0. 
By the chain rule, we rewrite the second term as JpPx, from which the group velocity 
u = J p. Here the subscripts t, x, p denote partial differentiation. The crucial feature 
is the inflection point in J(p), so that the group velocity can be either increasing 
or decreasing in p. We now employ the steady-state current ([2]) to determine the 
evolution of an initial density heterogeneity on length and time scales large compared 
to microscopic scales for the step initial condition 

p{x,t = Q) = r - (3) 

I p_i_ X > U . 

As sketched in Fig. HI the difference in the group velocity to the right and left of the step 
determines whether a continuous, discontinuous, or a composite density profile emerges. 
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It is worth noting that similar results for density profiles is obtained for an 
asymmetric exclusion process with another form of cooperative interactions |T8l [T9] . 
In these works, the same qualitative phase diagram as in Fig. HJis obtained, despite the 
rather different natures of the microscopic interactions in their model. This similarity in 
long-time behavior arises because our main results apply for any asymmetric exclusion 
process with sufficiently strong cooperative interactions, as indicated by an inflection 
point in J(p). 



Shock/Inverted Shock: A propagating shock wave arises whenever the group velocity 
on the left exceeds that on the right, u{p-) > u{p+). Qualitatively, the faster moving 
particles catch up to slower particles on the right and pile up in a shock wave, just 
as freely- moving cars suddenly slow down upon approaching a traffic jam. In the 
conventional ASEP, all upsteps evolve into a shock (S) wave. For the CE, in contrast, 
only upsteps where both initial densities are above the inflection point, p/ < p_ < p+, 
evolve into shocks (Fig. |5]). Here, exclusion is sufficiently strong that the group velocity 
is a decreasing function of density. Strikingly, a propagating shock wave also emerges 
from a downstep in CE when the initial densities are both below the inflection point, 
P/ > p_ > p+. In this regime, Jpp = Up > 0; that is, cooperativity is sufficiently 
strong that particles in the high-density region on the left have a greater group velocity 
and therefore pile up at the interface. We term this singularity an inverted shock (IS) 



For both shocks and inverted shocks, the density is given by the traveling wave 
proffie p = p(x — ct). We obtain the shock speed c by equating the net flux 
into a large region that includes the shock, J(p+) — J{p-), with the change in the 
number of particles, c(p+ — p_), in this region [22] to obtain the standard expression 
c = [J(p+) — J(p-)][p+ — p-]; this holds both for conventional and inverted shocks. 



P_ 




Figure 4. Phase diagram of tlie CE model for an initial density step {p-,p+), with 
Pj the inflection point in J{p). A typical density profile p{z) is sketched for each of the 
six regions: (R/IS) rarefaction/inverted shock, (R) continuous rarefaction, (S) shock, 
(C/S) compression/shock, (C) continuous compression, (IS) inverted shock. 



(Fig. ED. 
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Figure 5. (left) Evolution of an upstep for A = |: (C) continuous compression wave 
for p- = i, p+ = |; (C/S) composite compression/shock for p_ = |, p+ = (S) 
shock for /?_ = i, p+ (right) Evolution of a downstep for A = |: (R) continuous 

rarefaction for p_ = 1, p+ = (R/IS) composite rarefaction/inverted shock for 
p- = 1, p+ — |; (IS) inverted shock for p_ = 0.325, p+ = |. The dashed line is the 
locus Jp = z and the solid black curves are analytic predictions. Simulations are based 
on 10^ realizations up to i = 4 x 10^. 



Continuous Rarefaction/ Compression: A density step gradually smooths out when 
the when the group velocity to the left is less than that on the right, u{pJ) < u{p^). 
Here the faster particles on the right leave open space for the slower particles, similar 
to a cluster of stopped cars that slowly spreads out after a stoplight turns green. In 
ASEP, a downstep always evolves to a continuous rarefaction (R) wave. This continuous 
rarefaction also occurs in CE when both initial densities are above the inflection point, 
P- > p+ > Pi. At these high densities, exclusion dominates, as in the ASEP, which 
causes the group velocity to decrease with density. 

In striking contrast to the ASEP, an upstep can continuously smooth out in CE 
when the initial densities are below the inflection point, p_ < p+ < pj. In this 
regime, cooperativity is sufficiently strong that particles in the high density region on 
the right move faster than those on the left. Thus instead of a shock wave, a continuous 
compression (C) wave develops (Fig. [5]). We determine the density profile by assuming 
that it a function of the scaled variable z = x/t. Substituting p{x,t) = p{z) into the 
continuity equation gives —zp^ + JpPz = 0. Thus the profile consists either of constant- 
density segments (p^ = 0) or else z = Jp. Matching these solutions gives [HI [21] 

{p_ z < z_ , 

Hz) z.<z<z+, (4) 
p+ z> z+ , 

where I{z) is the inverse function oi z = J p. For a continuous profile, the cutoffs z_ and 
z+ are determined by matching the interior solution I{z) with the asymptotic solutions: 
I{z±) = p± or equivalently, z± = Jp{p±). 
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Composite Rarefaction/Compression and Shock: In CE, a continuous rarefaction 
or compression wave can coexist with a shock wave. This phenomenon occurs when 
the group velocity on the left is initially less than that on the right but also with 
the constraint that the initial densities lie on either side of the inflection point. 
Consequently one side of the step is in the exclusion-dominated regime and the other 
is in the cooperativity-dominated regime, or vice-versa. In particular, a composite 
rarefaction/inverted shock (R/IS) wave emerges from a downstep when > pj > p+, 
so that u{pS) < m(p_|_). As in the case of the continuous rarefaction wave, the downstep 
begins to smooth out from the rear. Consequently, cooperative interactions become more 
important as the density at the leading edge of this rarefaction decreases. Eventually 
this leading density reaches the point where the particle speed matches that at the 
bottom of the downstep and the rarefaction front terminates in an inverted shock. 

Correspondingly, an upstep can evolve to a compression wave with a leading shock 
when the densities satisfy p_ < p/ < p+ and u{pJ) < u{p+). In this case, the leading 
particles initially race ahead, leaving behind a profile where the density increases with 
X. However, this increase cannot be continuous because eventually a point is reached 
where the speed at the front of this continuous wave matches that of the top of the 
upstep. After this point, a pile-up occurs and a shock wave forms. We call this profile 
a composite compression/shock (C/S) wave (Fig. [5]). 

The functional forms of the composite rarefaction/inverted shock and composite 
compression/shock profiles are still given by Eq. (jl]), but the criteria to determine the 
cutoffs z± are now slightly more involved than for continuous profiles. The location of the 
left cutoff, is again determined by continuity, namely, = P- or, alternatively, 

= Jp{p^). To determine the right cutoff z+, note that in a small spatial region that 
includes the leading-edge discontinuity, the density profile is just that of a shock or 
inverted shock wave. Thus the equation for the shock speed is 

where g+ = I{z+) is the density just to the left of the discontinuity. (Note also that 
z^ = Jp{q+) by definition.) To justify ([5]), we use the conservation equation that the 
particle number in z+] equals the initial number plus the net flux into this region: 

J I{z)dz = -p-Z- + P+Z+ - J{p+) + J{p-) ■ (6) 

We recast this expression into ([5]), by making the variable change z = Jp{p) and using 
I{Jp{p)) = p to write the integral as /p+ p Jppdp, which can be performed by parts. The 
resulting expression readily simplifies to ([5]). 

In summary, a diversity of wave singularities arise in asymmetric exclusion with 
sufficiently strong cooperativity. The minimum requirement for these phenomena is 
an inflection point in the current-density relation J{p). This inflection point leads to a 
group velocity that is an increasing function of density for p < p/, a dependence opposite 
to that in the conventional ASEP. The resulting non-monotonic density dependence of 
the velocity causes an initial density upstep or downstep to evolve to: shock/inverted 
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shocks, continuous rarefaction/compression waves, or a composite profile with both 
continuous and discontinuous elements. 
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